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Abstract. We show central limit theorems (CLT) for the linear statistics of symmetric 
matrices with independent heavy tailed entries, including entries in the domain of attrac- 
tion of a-stable laws and entries with moments exploding with the dimension, as in the 
adjacency matrices of Erdos-Renyi graphs. For the second model, we also prove a central 
limit theorem of the moments of its empirical eigenvalues distribution. The limit laws are 
Gaussian, but unlike the case of standard Wigner matrices, the normalization is the one 
of the classical CLT for independent random variables. 

1. Introduction and statement of results 

Recall that a Wigner matrix is a symmetric random matrix A = {aij)ij=i^,,,^N such that 

1. the sub-diagonal entries of A are independent and identically distributed (i.i.d.), 

2. the random variables y/Naij are distributed according to a measure /i that does 
not depend on and have all moments finite. 

This model was introduced in 1956 by Wigner [57] who proved the convergence of the 
moments 
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when fi is centered with unit variance. Moments can be easily replaced by bounded con- 
tinuous functions in the above convergence and this convergence holds almost surely. As- 
sumption 2 can also be weakened to assume only that the second moment is finite. The 
fiuctuations around this limit or around the expectation were first studied by Jonsson [28] 
in the (slightly different) Wishart model, then by Sinai and Soshnikov [35] with p <^ N^^'^ 
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possibly going to infinity with A^. Since then, a long list of further-reaching results have 
been obtained: the central limit theorem was extended to so-called matrix models where 
the entries interact via a potential in |27], the set of test functions was extended and the 
assumptions on the entries of the Wigner matrices weakened in [71 El |29], a more gen- 
eral model of band matrices was considered in [2] (see also [291 E] fo^' general covariance 
matrices), unitary matrices where considered in [TH], and Chatterjee developed a general 
approach to these questions in [T7], under the condition that the law /x can be written as 
a transport of the Gaussian law. Finally, but this is not really our concern here, the fluc- 
tuations of the trace of words in several random matrices were studied in [161 ESI [321 [26] . 
It turns out that in these cases 

Tt{AP) -E[Ti{AP)] 

converges towards a Gaussian variable whose covariance depends on the first four moments 
of /i. Moments can also be replaced by regular enough functions and Assumption 2 can 
be weakened to assume that the fourth moment only is finite. The latter condition is 
however necessary as the covariance for the limiting Gaussian depends on it. The absence 
of normalization by \/N shows that the eigenvalues of A fluctuate very little, as precisely 
studied by Erdos, Schlein, Yau, Tao, Vu and their co-authors, who analyzed their rigidity 
in e.g. [2111221 Eg. 



In this article, we extend these results for two variations of the Wigner matrix model 
where Assumption 2 is removed or weakened: we consider cases where /i does not have any 
second moment or cases where fi depends on A^, with moments growing with A^. In these 
cases, Wigner's convergence theorem (JT]) does not hold, even when moments are replaced 
by smooth bounded functions. The analogue of the convergence was studied when the 
common law jj, of the entries of y/N A belongs to the domain of attraction of an a-stable 
law (in which case the normalization in \/N has to be adapted) or /i depends on A^ and 
has moments blowing up with A^. We describe these models below more precisely. 

Definition 1.1 (Models of symmetric heavy tailed matrices with i.i.d. sub-diagonal en- 
tries). 

Let A = {aij)ij=i^,„^N be a random symmetric matrix with i.i.d. sub-diagonal entries. 

1. We say that A is a Levy matrix of parameter a m ]0, 2[ when A = X/a^ where the 
absolute values of the entries Xij of X are in the domain of attraction of a-stable 
distribution, more precisely 

(2) P{\X,,\>u) = ^ 

ill 

with a slowly varying function L, and 

aN = inf{M : P i\Xij\ > u) < ^} 
(a^ = L{N)N^''^, with Z(-) a slowly varying function). 
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2. We say that A is a Wigner matrix with exploding moments with parameter 
(C*fc)fc>i whenever the entries of A are centered, and for any k >1 



with Y.k>2 Cfe > 0- 

Notice that if there exists a > for a k > 2, then actually Ck > for any k > 2 
[3T| Appendix A]. A Levy matrix whose entries are truncated in an appropriate way is a 
Wigner matrix with exploding moments [lOl |3T1 [38] • The adjacency matrix of an Erdos- 
Renyi graph, i.e. a matrix A such that 

(4) Aij = 1 with probability p/N and with probability 1 — p/N, 

is also an exploding moments Wigner matrix. In this case the fluctuations were already 
studied in [M]. 

The weak convergence of the empirical eigenvalues distribution of a Levy matrix has 
been established in [TU] (see also [151 El US]) where it was shown that for any bounded 
continuous function /, 



N-^oo Jy J 

where fia is a heavy tailed probability measure which depends only on a. Moreover, /i^ 
converges towards the semicircle law as a goes to 2. 

The convergence in moments, in expectation and in probability, of the empirical eigen- 
values distribution of a Wigner matrix with exploding moments has been established by 
Zakharevich in |3B]- If "we assume that the sequence {Ck)k>i has at most exponential 
growth, i.e. that for a constant C, for all k, 



where /xc is a probability measure which depends only on the sequence C := {Ck)k>i- 

We shall first establish the fluctuations of moments of a Wigner matrix with exploding 
moments around their limit, namely prove the following theorem. 

Theorem 1.2 (The CLT for moments of Wigner matrices with exploding moments). 
LetA=( ^i,j)i,j='^,--;N be a Wigner matrix with exploding moments with parameter {Ck)k>i- 
Then the process 



(3) 





a.s. 




a.s. 



(6) 




converges in distribution to a centered Gaussian process. 
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The expression of the covariance of the process is of combinatorial nature, given in 
Section O Formula (133|) and Theorem 12.21 

Note that the speed of the central limit theorem is N'^/"^ as for independent integrable 
random variables, but differently from what happens for standard Wigner's matrices. This 
phenomenon has already been observed for adjacency matrix of random graphs [HI ED] 
and we will see below that it also holds for Levy matrices. It suggests that the repulsive 
interactions exhibited by the eigenvalues of most models of random matrices with not too 
large entries (see the first paragraph of this introduction) are absent or much less strong 
in the context of heavy tailed matrices. 

For Levy matrices, moments do not make sense and one should consider smooth bounded 
test functions. We start as is common in random matrix theory, with the study of the 
normalized trace of the resolvant of A, given for z G C\M by 

By the previous results, for both Levy and exploding moments Wigner matrices satisfying 
(|5]), A^~^ Tr G{z) converges in probability to a deterministic limit as the parameter tends 
to infinity. We study the associated fluctuations. 

In fact, even in the case of Wigner matrices with exploding moments, the CLT for 
moments does not imply a priori the CLT for Stieltjes functions even though concentration 
inequalities hold on the right scale, see [H]. Indeed, one cannot approximate smooth 
bounded functions by polynomials for the total variation norm unless one can restrict 
oneself to compact subsets, a point which is not clear in this heavy tail setting. However, 
with additional arguments based on martingale technology, we shall prove the following 
result, valid for both Levy matrices and Wigner matrices with exploding moments. 

Theorem 1.3 (CLT for Stieltjes transforms). Let A be an N by N Hermitian matrix with 
i.i.d. real entries below the diagonal defined in Definition \l.l\ satisfying ([5]) in the case of 
Wigner matrices with exploding moments. Let 

(7) 0Ar(A) :=E[exp(-zA|aii|2)], 

defined for any A in C^, with := {A G C ; Q'A < 0}. Assume that there is a function $ 
on C~ such that as N ^ oo, we have the convergence 

(8) iV(0iv(A) - 1) ^ $ (A). 
uniformly on compact subsets of . Then, the process 

(9) ( ^TtG(z) - ^E\TtGU 

converges in distribution (in the sense of finite marginals) to a centered complex Gaussian 
process Z{z), z G C\R, whose covariance Cov (^Z (z) , Z {z')) := K[Z{z)Z{z')~\ = G{z,z') is 
given in Formulas (I13p (1141) and (llSp below. 
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Note that the uniform convergence on compact subsets, in ([8j), imphes that $ is analytic 
on C~ and continuous on C^. 

As a consequence, we can extend the central limit theorem to a larger class of functions. 

Corollary 1.4. Let A be the set of analytic functions f on the strip {z = a + ib,a & 
— ^} f^'^ some e > so that there exists 7 > 1 so that for a large enough 

(10) \f{a±ie)\<a-^ . 
Then under the hypotheses of Theorem \l.'J\ the process 

(11) (^^(/) := ;^Tr/(A) - -^E[Tr/(A)] ■ f ^ A) 

converges in law towards a centered Gaussian process {Z{f) ; f ^ A) with covariance de- 
fined by the facts that Z{f ) = Z{f), with f the function of A defined by the formula 
Tiz) ■= fiz), and 

Cov {Z{f),Z{g)) = -1- [dz [ dz'f{z)g{z')C{z,z'), 

where T =] — 00 + ze, +00 + ze[ U ] + 00 — ze, —00 — ie[ and C{z, z') is the function of 
Theorem \1.3[ The result extends to the closure of A for the total variation norm. 

Examples: before giving the explicit form of the covariance, let us first discuss the hy- 
pothesis made at Equation (|8]). 

• For a Levy matrix, ii A = X/aj^ with \Xij\ in the domain of attraction of an a- 
stable law, Xfj is in the domain of attraction of an a/2-stable law, and we have 
(see [23]) that holds with 

(12) $(A) = -(t(zA)"/2 

for some constant cr > (in this text, as specified in the notations paragraph at 
the end of this section, the power functions have a cut on R~). 

• For a Wigner matrix with exploding moments, at least on the formal level, we have 

N{M>^)-l)=Y,t^NE[\ann. 

k>l 

Therefore, if (E]) holds together with the reinforcement of into 

N>1 

then OS} is satisfied with $(iA) the exponential power series of the C^'s, i.e. $(A) = 

^ (_iA)^ 

• In particular, since the adjacency matrix A of an Erdos-Renyi graph defined by (jl]) 
satisfies (E]) with Ck = p for any k > 1, then ([S]) holds with 

$(A) =p{e-'^-l). 
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The function C{z, z') in Theorem 11.31 is given by C{z, z') = L{z, z') — L{z)L{z'), where 
(13) L{z) = / -dJ'^+P^^'^dt, 



t 

rl /•sgnS{^)oo /"Sgn 3(2)00 -r 

The maps pz and p"^, are given by p"2/(t, t') = y [$(ta + t'a')] , p2(t) = p"^/(t, 0), where 
(a, a') is distributed according to the non random weak hmit asA^— )-oo, k/N ^ u oi 



2 ..*.+.*'.'+p:,,,(M')d^d^/d^_ 



where Gk{z) = [z — A^)"^ for the matrix obtained from A by deleting the k-th row and 
column, and G'j^{z') = {z' — A'i^)~^ for A'/^ a copy of Ak where the entries (i, j) for i,j>k are 
independent of Ak and the others are those of Ak- Also, we assume tQz > and t"^z' > 
in (1151). 



The existence of the limit (1T5|) is a consequence of a generalized convergence in moments, 
namely the convergence in distribution of traffics, of {Ak^A'/.) stated in [31], see Lemma 
13.21 However, under stronger assumptions, an independent proof of this convergence and 
an intrinsic characterization of p" ^, are provided in Theorem 11.61 below. 

Let us first mention that the map pz is the almost sure point-wise limit 
(16) Pzit) = 

and it characterizes the limiting eigenvalues distribution of A. Indeed, under the assump- 
tions of Theorem 11.31 for all z G C"*", we have the almost sure convergence 



oo 



(17) lim — TrGAr(^) = i / e'^'^f^^'Mt. 

N^oc N Jq 

This fact was known in the Levy case [10] and is proved in Corollary 13. 31 for Wigner matrices 
with exploding moments. Let us first give a characterization of pz- 

Theorem 1.5 (A fixed point equation for Pz{t))- We assume the hypotheses of Theorem 
\1.3\ and, moreover, we assume that we have the following decomposition for all z G C~ .■ 



(18) $(z) = / g{y)e'Uy 

Jo 

where g{y) is a function such that for some constants K,'y > —1, n > 0, we have 

(19) \g{y)\ < Kly<iy^ + Kly>iy^, > 0. 



Then the function {z G C~^, A G M"^) H- Pzi^) is analytic in its first argument and contin- 
uous in the second, with non positive real part, and characterized among the set of such 
functions by the formula 

/■oo 

(20) p,(A) = A / g{\y)e'y^^''^^y^dy. 



Examples: Assumption (|T8i) is satisfied for our main examples. 

• In the heavy tailed case, we have $(a;) = —a{ix)^, hence Hypothesis f fTSj) is true, 
with g{y) = Cy^~^. Indeed, for any z G C"^, an application of residues formula 
gives 

(21) r(a/2) = -i {-izt)'^-^e'''zdt. 

Jt=o 

• In the case of the Erdos-Renyi adjacency matrix, $(x) = p{e^^^ — 1) so that 

POO _^ 2 r°° 1 

(22) <^(x) = -p / Ji(s)e'^ds = -p / Ji(2v^)e""'^^ds 

Jo Jo Vs 

where Ji{s) = | ^fc>o ^k\^k+iy. ^ Bessel function of the first kind. Ji is bounded 
on M+. 

• More generally, if $ is so that $(x~^) is in L^, by Fourier integral formula we can 
write 



with $_i the Fourier transform of $(l/x). Hence, if $_i is supported on we 
have 

POO 

^x) = / <l>_i(i/)e-^^/^dy, 
Jo 

so that our hypothesis is satisfied if $_! satisfies (HM . 

Note that the asymptotics of Wigner matrices with bounded moments is also described 
by (12U]) . In this case $(A) = — iA, so g{y) = —i and Pz(t) = —it lim jjTiG^z), which 

N^oo^^ 

leads, by formula (|T7|) . to the classical quadratic equation 
(23) s{z) := lim IzTtG{z) ^ 



N^ooN Z - s{z) 

Let us now give a fixed point characterization for the function Pz,z'{t,t') of (fT5i) . 

Theorem 1.6 (A fixed point system of equations for ^^lit, t')). We assume the hypotheses 
of Theorem \l.'J\ and we assume one of these assumptions. 
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1) $ satisfies the formula 



(24) $(x + ?/) = j j e'-^' y dT{v,v') + j e'~dfi{v) + j e' y dfi{v') 

for some complex measures r, fi on respectively x and such that for all 
b > 0, j e~'''"d\fi\{v) is finite and for some constants K > 0, —1 < < and k> 0, 
and 

(25) "^'^[^^^^ ^ < J^(t;n„e]o,i] + t;"1.6]i,oo[) (^/^l.'6]o,i] + t^'"l.'e]i,oo[) • 
2) A is a Levy matrix, so that $(x) = — (t(zx)°/^. 

Then, the parameter p^^z' of f lT^ zs g'zwen &y 

= ^Pri(^' + (1 - '")p";f'(^, s) 

where the pair {p^'li{t, s), p^''^,{t, s)) zs the unique solution which are analytic on A = {z, z' : 
tQz > 0,s'^z' > 0} and uniformly continuous on compacts in the variables {t,s) (and /3- 
H older for (5 > a/2 in the Levy matrices case, see Lemma \5. 1\) of the following fixed point 
system of equations 

fsgn, CO /"Sgn./ oo , , 

p:'l,{t,s) =11 e^T^+^T^'e<.'(T'^)dr(^;,t;') 



^0 

sgn. oo fsgn^i oo 

e'hf~~^^-Mp{v) + / e'^-eP^'^^-^dp{v) 
Jo 

Jo Jo 

rsgn^ oo /"Sgn^/ oo 

(26) + / e'h''^^^-Mp{v) + / e'^-eP^'^^-Mp{v), 

Jo Jo 

with the notations sgn^ := sign(532), sgn^ := sign(532;'). 

Examples: in our main examples: 

• is satisfied in the case of the Erdos-Renyi graph adjacency matrices. Indeed, 
in this case, by (1221) . 

$(x + ?/) = p(e^+s^-l) =p(e^'-l)(e^-l)+p(e^-l)+p(e^-l) 
= p Ji{v)dv / Ji{v')dv'e^^~^^^ 



l2 



—p / Ji{v)dve^*^ —p Ji{v)dve 
Jo Jo 

so that if u denotes the measure on defined by 

du{v) = Ji(2v^)^dw 
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then we can take r = pz/®^ and fi = —piy, which satisfies our assumption with 
7 = -1/2, K = 0. 

• For Levy matrices, where ^{x) = CaX°'^'^ with Ca = — cz"/^, fl2^ holds as welL We 
shall use later the fact that in this case, for all x,y E (with a constant Cq that 
can change at every line), 

(27) $(x- + y-^) = C^{\ + i)"/^ = C^j^^j^^{x + yT'- 

POO POO POO 

= C„ / dw dw' / di;y;"/2-l^^/)a/2-l^-a/2-lg™xW2/^g»(x+j/) _ ^) 

io io Jo 

(where we used the formula 2;°^^ = Ca fa/2+i ^^■^ ^ ^ ^ 

(0,2), which can be proved as (1211) with the residues formula) so that holds 
with /i = and t{v,v') with density with respect to Lebesgue measure given by 

r+oo 

(28) Ca / - ur/'-\v' - ur/^-'lo<u<vAv' - v''/^~\vr^^~'}du. 







Unfortunately r does not satisfy (125!) as its density blows up at f = v': we shall 
therefore prove the uniqueness of the fixed point system of equations separately. 



Let us conclude this introduction with three remarks. 

(1) Let A = X/ttN be a Levy matrices defined at Definition II. ip but with a = 2. Then 
using Example c) p. 44. of [23] instead of the hypothesis made at Equation ([7j), 
one can prove that as — 00, the spectral measure of A converges almost surely 
to the semi-circle law with support [—2, 2] (see fl2^ ). This result somehow "fills the 
gap" between heavy-tailed matrices and finite second moment Wigner matrices. It 
allows for example to state that if P(|Xjj| > u) ~ cu~^, with c > 0, even though the 
entries of X do not have any second moment, we have that the empirical spectral 
law of , ^ converges almost surely to the semi-circle law with support [—2, 2]. 

yJcn\og{n) 

(2) About applications of our results to standard Wigner matrices, i.e. real sym- 
metric matrices A with i.i.d. centered entries (up to the symmetry) with variance 
1/N : one can see from the proofs that Theorems 11.31 and 11.51 and Corollary 11.41 
also apply for such matrices and that the function $ of ([8]) is then linear, which 
implies that L{z,z') = L{z)L{z') for all 2,-2', so that the covariance is null. In 
this case, (!20|) is the self-consistent equation satisfied by the Stieltjes transform of 
the semi-circle law, namely (1231) . and Corollary 11.41 only means that for functions 
/ G v4, we have, for the convergence in probability. 



(29) 



Tr/(A)-E[Tr/(A)] =o(v^). 
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When the entries of the Wigner matrix A have a fourth moment, (l29i) does not 
bring any new information, because it is aheady well-known that its LHT converges 
in distribution to a Gaussian law. However, for Wigner matrices whose entries 
have a second but not a fourth moment, (129|) brings new information. Indeed, 
for such matrices, which could be called "semz heavy-tailed random matrices" , the 
convergence to the semi circle law holds (see [1] or the remark right above that 
one) but the largest eigenvalues do not tend to the upper-bound of the support 
of the semi-circle law, are asymptotically in the scale (with a G (2, 4) as in 

Equation ([2]) when such an exponent exists) and distributed according to a Poisson 
process (see [3]), and it is not clear what the rate of convergence to the semi-circle 
law will be. Equation fl2^ says that this rate is ^ N'^/"^ . 
(3) About recentering with respect to the limit instead of the expectation, it depends 

on the rate of convergence in ([3]) or in ([7]) . For instance, if N'E,{a'^)—Ck = o{-\fN ) 
for any > 1, then 





— TiA^ 


- lim E 


— TiA^ 




IN J 




IN J 



0, 



but otherwise a non trivial recentering should occur. See the end of Section [2l 



Notation: In this article, the power functions are defined on C\M_ via the standard 
determination of the argument on this set taking values in (— vr, tt). The set C"*" (resp. C~) 
denotes the open upper (resp. lower) half plane and for any z G C, sgn^ := sign(53z). 



2. CLT FOR THE MOMENTS OF WiGNER MATRICES WITH EXPLODING MOMENTS 

The goal of this section is to prove Theorem 11.21 In order to prove the CLT for the 
moments of the empirical eigenvalues distribution of A, we use a modification of the method 
of moments inspired by [SU] which consists in studying more general functionals of the 
entries of the matrix (the so-called injective moments) than only its moments. We describe 
this approach below. 

Let A be a Wigner matrix with exploding moments. Let > 1 be an integer. The 
normalized trace of the i^-th power of A can be expanded in the following way. 

1 ^ 

— A{ii,i2)...A{iK-i,iK)A{iK,H) 

ii,...,iK=l 

^ ^ ^ A{ii, za) . . . A{iK-i, iK)A{iK, ii), 

^ V ' 

where V{K) is the set of partitions of {!,... ,K} and S.,^ is the set of multi-indices i = 
(?!,..., ix) in {1, ... , A^}-^ such that n m ^ i^ = im- 



N 



TrA 



K 
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We interpret as a functional on graphs instead of partitions. Let vr be a partition of 
{1, . . . ,K}. Let T'" = (y,E) be the undirected graph (with possibly multiple edges and 
loops) whose set of vertices V is vr and with multi-set of edges E given by: there is one 
edge between two blocks Vi and Vj of vr for each n in {1, ... , K} such that n & Vi and 
n + 1 & Vj (with notation modulo K). Then, one has 

(30) r° (tt) = r° n 
if for graph T = (V, E), we have denoted 

(31) ^^^^)=^ E n^('^(^))' 

<p:V^[N] ee-E 
injective 

where [A^] = {1, . . . , N} and for any edge e = {i,j} we have denoted A (0(e)) = A(^(f){i), 0(j)) . 
There is no ambiguity in the previous definition since the matrix A is symmetric. 

In order to prove the convergence of 



(^^(^0),,>i=(^TrA^-E 



K>1 



■k<^VJkV{K) 



Tr 

to a Gaussian process, it is sufficient to prove the convergence of 

(32) {ZN{T-)).,.,r,K) ■■= [^{rUn " E[r^[Tl])) 
to a Gaussian process, since 

(33) {Zr,{K))^^^=( J2 ZNiT^) • 

Before giving the proof of this fact, we recall a result from [31], namely the convergence of 
r^[T'^] for any partition vr. These limits are involved in our computation of the covariance 
of the limiting process of (Zjv(T''))^^^^^^^^, and this convergence will be useful in the 
proof of the CLT for Stieltjes transforms latter. 

Proposition 2.1 (Convergence of generalized moments). Let A he a Wigner matrix with 
exploding moments with parameter {Ck)k>i- For any partition n in UkV{K), with t^[T'^] 
defined in (15U]1 . 

^ Ar_s.oo \ otherwise, 

where a fat tree is a graph that becomes a tree when the multiplicity of the edges is forgotten, 
and for T such a graph we have denoted qk the number of edges ofT with multiplicity 2k. 

Theorem 2.2 (Fluctuations of generalized moments). Let A be a Wigner matrix with 
exploding moments. Then, the process (^Zj^fiT'^)^^^^ v(k) ^^/^'^^^ (133) converges to a 
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tered Gaussian process 



centerea uausstan "process \ z(i " ) ) , whose covariance is given by: for any T'^^, T'"^ , 

where t^[T] is given by Proposition \2.1\ and Pn(T'^i, T'^^) is the set of graphs obtained by 
considering disjoint copies of the graphs T'^^ and T'"^ and gluing them by requiring that 
they have at least one edge (and therefore two "adjacent" vertices) in common. 



Proof. We show the convergence of joint moments of (^Zn(T'^)Y Gaussian distribution be- 
ing characterized by its moments, this will prove the theorem. Let Ti = {Vi, Ei), . . . ,Tp = 
{Vn, En) be finite undirected graphs, each of them being of the form T'" for a partition vr. 
We first write 



E 



Zn{Ti) . . . ZN{Tn) 



1 



N2 



E 

'1, .■■></>« 

4>i:Vj^[N\ inj. 



E 



eeE, 



V 

a;jv(</'iv,0n) 



E 



N2 
Vn) (</-! 



where 



V{Vi, . . . , Vn) is the set of partitions of the disjoint union of Vi, 
contain at most one element of each Vj, 

So- is the set of families of injective maps, (pj : Vj — )■ [N],j = 1, 
any v G Vj, v' G Vj', one has = (t>j'{v') v v'. 



Vn whose blocks 
, n, such that for 



First, it should be noticed that by invariance in law of A by conjugacy by permutation 
matrices, for any a in V{Vi, . . . , Vn) and . . . , (pn) in Sa, the quantity ujN{4>iy • • • ; 4>n) 
depends only on a. We then denote ujn{o') = (^n{(Pi, ■ ■ ■ , (pn)- Moreover, choosing a 
partition in V{Vi, . . . , Vn) is equivalent to merge certain vertices of different graphs among 
Ti, . . . ,Tn. We equip ViVi, . . . ,Vn) with the edges of Ti, . . . , T„ and say that two vertices 
are adjacent if there is an edge between them. We denote by V^{Vi, . . . ,Vn) the subset 
of V{Vi, . . . , Vn) such that any graph has two adjacent vertices that are merged to two 
adjacent vertices of an other graph. By the independence of the entries of X and the 
centering of the components in un, for any a in ViVi, . . . , Vn) \ V^iVi, . . . , Vn) one has 
|^Ar(o") = 0. Hence, since the cardinal of S„ is jj^^^, we get 



E 



(Ti)...Z;v(T„ 



{N-\a\)\ 



(34) 



<7ePj(yi,...,y„) 



a£Vi(Vu...,Vn) 
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Let o in Vf^iVi, . . . , Vn). We now analyze the term un{o-). We first expand its product. 

W7v(a) = 



(-l)"-l^lE 

BG{l,...,n} 



jeBeeEj ^ j^B '-eG-Ej 



Let B C {1, . . . ,n}. Denote by Tb the graph obtained by merging the vertices of Tj, 
j E B that belong to a same block of a. For any k > 1 denote by pk the number of vertices 
of Tb where k loops are attached. For any k > £ > 0, denote by qk/ the number of pair 
of vertices that are linked by k edges in one direction and i edges in the other. Denote by 
/iTV the common law of the entries of \/N A. By independence of the entries of A, for any 
(01, . . . , (pn) in So-, one has 



E 



n n Mu^)) 

'-jeBeeEj 



n 

k>l 



N2 



(35) 



N-\Eb\ -Q 



k>l 



n 

k,e>o 



I\ 2 



n 

kl>0 



/t^+M/i^(t) 



where \Eb\ is the number of edges of Tb once the multiplicity and the orientation of edges 
are forgotten. Recall assumption for any k > 1, 



NE[{aijf'] =E 



By the Cauchy-Schwarz inequality, for any A; > 1, 



tMflNit) = 0(1). 



Hence, since the measure /Uat is centered, the quantity Sn{B) is bounded. Denote T^- the 
graph obtained by merging the vertices of Ti, . . . , T„ that belong to a same block of a, \E(j\ 
its number of edges when orientation and multiplicity is forgotten, and by Co- its number 
of components. We obtain from (134|) and (135|) 

E 



(Ti)...Zjv(T„ 

J2 Yl iV-^+l'^l-l^^l-^^^^ 1%}! (-l)"-l^l<5^(S) J] 5Ar({j}) (1 + 0(iV-^)) 

(T67'tt(V^i,-,V„)Bc{l,...,n} i^B 

(TePB(yi,...,y„)Bc{i,...,n} 

x(_i)"HBl5^(5) J]5^({j}) X (l + 0(iV-i)) 
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A partition a G V^iyx, . . . , Vn) induces a partition (t of {1, ... , n\: i j if and only if 
Ti and Tj belong to a same connected component of T^-. Denote by V2{n) the set of pair 
partitions of {1, . . . , n}. One has 

(36) iV^''-t = W,(„) + 0(iV"i). 

Secondly, one has l-Eo-l — \Eb\ — J2ji^B 1-^0} I — ^ with equality if and only if i? = 
{1, . . . ,n}, so that 

(37) n\^^\-\Eb\-J:,<^b I%}I = l^^ii^...^,} + 0{N-^). 

Moreover, by the relation [211 Lemma 1.1] between the number of vertices, edges and 
components of a graph, \a\ — — \E„\ is the number of cycles of T^, the graph obtained 
from Tfj by forgetting the multiplicity and the orientation of its edges. Hence, we have 

(38) ivH-c.-|^.| = %i3,f_,, + 0(iV-^). 

By (13SD, dSZ]) and ([3SD, if we denote 6n{(t) = 5jv({l, ...,n}) we get 
E 



Zn{Ti) ■ ■ ■ Z^iTn 

If,, is a forest of 2 trees 5Ar(cr)+0(A^ ^) 



7re7'2(n) <j6Pj(yi,...,y„) 
s.t. t7=7r 



(39) = E n E %isatree ^7v(^)+0(iV-l), 

where we have used the independence of the entries of A to split 6^. The case n = 2 gives 



E 



aeVt{Vi,V2) 



M(2)(Ti,T2) 

where (5(cr) = lim 5n{c), which exists since 7^ is a tree. Indeed, in the definition fl35i) of 

A''— 5>oo 

5(cr), we have pk = Qk/ = for any k ^ i. Moreover, we obtain that 6{a) = t^[T„] defined 
in Proposition l2.lt and the sum over a on Vf^iVi, V2) can be replaced by a sum over graphs 
T obtained by identifying certain adjacent vertices of Ti with adjacent edges T2, since r°[T] 
is zero if T is not a fat tree. We then obtain as expected the limiting covariance 

M(2)(Ti,T2)= Yl ^°[^]- 

TePB(Ti,T2) 



The general case n > 3 in fl39l) gives the Wick formula 

E[z^{T,)...Z^{T^)j= Y II M^^\Ti,T,) + oil), 

7rG'P2(n) {jj'}67r 

which characterizes the Gaussian distribution. □ 
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Remark about the recentering with respect to the hmit instead of the expectation: in 
our computations, up to (139|) the errors terms are of order 0{N~^), and so if A^E[a^^] = 

Ck + o{\/N ) for any k > 1, then 5iv(o") = 6{a) +o{\/N ) and hence one can replace the 

0(1) by o{\/N ) at the end of the proof. Similarly, we would find that 



E 


— TiA'' 


= hm E 






IN \ 




IN J 



so that we would obtain the same CLT if we recenter with the limit or the expectation, as 
noticed in the introduction. 



3. CLT FOR Stieltjes transform and the method of martingales 

Let A = {aij)ij=i^,,,^N be either a Levy matrix of parameter a G (0,2), or a Wigner 
matrix with exploding moments of parameter {Ck)k>i- For any 2 in C \ M, recall that 
G{z) = [z — A)~^. To prove Theorem 11.3^ we show that any linear combination of the 
random variables 



1 



Zn{z) := ^TtG{z) -E -^TtG{z) , zeC\ 



1 



N 

and their complex transposes converges in distribution to a complex Gaussian variable. 
Since G{z) = G{z), it is enough to fix a linear combination 

p 

(40) M{N) ■.= J2^^Z,,{zi), 

1=1 

for some fixed p > I and Ai,...,Ap G C, zi,...,Zp G C\]R, and prove that M{N) is 
asymptotically Gaussian with the adequate covariance. 

3.L The method of martingales differences and reduction to the case p = 1. For 

all A^, we have 

AT 

M{N) - E[M{N)] = J2 Xk with Xk := (Efc - E,_i) [M(iV)] , 

k=i 

where E^ denotes the conditional expectation with respect to the cr-algebra generated by 
the k first columns of A. In view of Theorem 16.71 of the appendix about fluctuations of 
martingales, it is enough to prove that we have 

N p 

(41) j]A.A,C(z.,z,), 

k=l j,i=i 

N p 

(42) Y.^k-i[\X,\'] 5^A.A;C(^„z-) 
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and that for each e > 0, 



(43) Y.^[\X,\H\x,i>s] 0. 

k=l 

Notice first that (HTj) imphes (1421) . Let us now prove fl43l) . The proof of fl4ip will then 
be the main difficulty of the proof of Theorem 11.31 

Proof of (143|) . Let be the symmetric matrix with size A^ — 1 obtained by removing the k- 
th row and the k-th column of A and set Gk{z) := jza^- Note that 'K^Gkiz) = Kk„iGk{z), 
so that we can write 

Xk = 5^ A, X (E, - Efc_i)[^(Tr - Tr 



1=1 * 
Hence by (194|) of Lemma [6.51 in the appendix, there is C such that for all and all k, 

Cmaxj |Ai||S2i|"^ 



'N 

Thus for N large enough, we have that for all /c, |^fcpl|Xt.|>e = and fH51) is proved. □ 



We now pass to the main part of the proof of Theorem II. 3^ namely the proofs of f HT]) 
and f l42p . It is divided into several steps. 



Let us first show that we can get rid of the linear combination in fl40p and assume p = 1. 
We have 

P 1 

Xu = (Efc - Efc_i) [M{N)\ = J2 A.(Efc - Efc_i) Tr G{z^)], 



hence, as Ti G{zi) = Ti G(zi), both X| and are linear combinations of terms of the 

form 

(Efe-Efc_i)[4=TrG(^)] X (E^ -Ej,^,)[^TtG{z% 



with z' G C\M. As a consequence, we shall only fix 2;, z' G C\M and prove that for 

n := (Efc -E,_i)[4=TrG'(z)] and Y^' := (E,, - E,_i)[^ Tr G'(z')], 
V A/ V A* 



we have the convergence in probability 

TV 

(44) G^ := ^E,.i[nF,:] G{z,z'). 



k=\ 
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First, for Gk as introduced in the proof of (143|) above, as ¥.kGk{z) = Kk-iGk{z) again, 
we have = (E^ — Kk-i)[^^{TT G{z) — TiGkiz))]. Hence by Lemma [63| we can write 



(45) 



1 l + a*Gfc(2)2afc 



N 



z- ttkk- alGk(z)ak ' 



where is the k-th column of A where the A;-th entry has been removed. To prove (jSj), 
we shall first show that we can get rid of the off diagonal terms '^jj!,i^k{j)^k{^)Gk{z)j£ 
and ^j^£Sik{j)sik{i){Gk{z)'^)ji in the above expression. 

3.2. Removing the off-diagonal terms. In this section, we prove that we can replace 
Yk in dSD by 



i + E7=r>,(jr(G,(z) 
^ -j2f=i^kUyGk{z)_i_i 



(46) 



Y, :-- 



N 



- Ek-i] 



Note first that by Equation fl97j) of Lemma [6.61 in the appendix, we have, with probability 
one, for all z G C\]R, 

'i + E.^kUnOkizYh 



z - T^i^kUyOkiz] 



< 2 



hence 
(47) 



o^i-i 



\Y,\ < 



N 



Lemma 3.1. Let Y^ be defined in the same way as Yk in fH6l) . replacing z by z' . Set also 



N 



(4J 



CN:=Y,^k-i[Ykn]. 



k=l 



Then for Gn defined as in f l44p . as N ^ oo, we have the convergence 



(49) 



Cn — Cn 



0. 



Proof. First note that for both Levy and exploding moments matrices, for any (A^ — 1) x 
(A^ — 1) matrix B with spectral radius 



where there exists 5 > and a finite constant G so that for all t > 0, all matrix B, 



(50) 



P(k7v(fi)| >t) < 
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By definition, Sn^B) = ^k(y^)^k{j)Bjk. In the case of exploding moments, we liave 
E |£7v(5)|'] = ^^E[afc(z)afc(j)%afc(z')afc(/fe] 

< 2 5^E[a,,(^)2]E[a,,(jr]|6, 



< ^(iVE[a,(lf])^||fi|'^ 



from which flSUl) follows by Tchebychev's inequality with 5 = 1/2. For Levy matrices, we 

put for K G (0, 2(23^) 

£n{B) = ^ afc(£)afc(j)_B^j = (Imaxj |afe(j)|<Ar« + Imaxj |afc(j)| >Ar'' )^Af • 

By [TUl Lemma 4.3], and the same computation as above, for any t > we have 

2 



P ({k^l > n {max \a,,\ < iV^) < 



(2-a)C-^2 



whereas there exists a finite constant C so that 



p(maxK,|>iV'^) <^ 



which prove flSU]) in this setting with 26 = min{ttK, 1 — 2k(2 — a)}. 

We apply f lSUjl with i? = Gk{z) and -B = Gkiz)"^. Since the function of two complex vari- 



ables ip{xi, X2) := has a uniformly bounded gradient on the set {(xi, 0:2) G ; '^X2 < 
0}, it follows that 



1 + a*Gfc(^)2afc 1 + afc(jf ^ 



N 



where \e'j^\ < G'\\eN{Gk{z)) \ + \eN{Gk{zY)\] for some finite constant G' . Now, we can use 
dnS]), (Ell) ^"^^ dSZl) to assert that the two first terms, in (EI]), are bounded by a constant 
uniformly in N and k. Hence, e'j^ is uniformly bounded. It follows that there exists a finite 
constant G" so that for any t > 0, 



Taking t = N ''we conclude 

E[i4r] < 



, ,2, ^ C""logiV 



To sum up, we have 
(52) n = n + 7. n = n + 7fc, 

where the error terms Jk^li 0{N^^~^ x logA^) uniformly in k for the norm. 
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But 

A, := Y,Yl - Y,Yl = Y^^ + ^,'7^ + 7,7^, 

so that 

E[|Afe|] < \\%\\lA\i'Al- + ||n'llL2||7fc||L2 + hfclUHlTfcllL^, 
with Ffc = Oloo(A^-i/2) (and the same for Y^ by (14 7p so that (152 p imphes that 

which proves the lemma. □ 

It remains to prove that the sequence introduced at (HHI) converges in probabihty as 
goes to infinity. Note that we have 

(53) C^=l iVE,_i[n(„)l>ydu with A;(n) := [iVn]. 

By (H7|) . the integrand is uniformly bounded by 4|Q'2;|~^. Hence, by dominated convergence, 
it is enough to prove that for any fixed u G (0, 1), as /c, tend to infinity in such a way 
that kjN — J- M, we have the convergence in probability of 

(54) iVEfc.iiny,']. 
Now, set 

^ l + E,a.(jr(G.(^)^),, ^ l + E,a.(jr(G,(^0%- 

^ ^ ^-E,a.(jm(.),, ^'^^ .-Ef=ia.(jTG,(z'),/ 

We have ^J~NY\, = (E^ — Efc_i)/fc and the same for Y^,, so 

iVEfe„i(ffey;') = E,_i[Efc/fcEfc/^] - Efc_i/A_i/^. 

Let us denote by Ea^. the expectation with respect to the randomness of the fc-th column 
of A [i.e. the conditional expectation with respect to the a-algebra generated by the aj^'s 
such that k ^ {i,j}). Note that ^k-i = °'^k = ° Ea^,, hence 

iVEfc_i(ny;0 = EajE,/fcE,/;] - EfcEaJfc X EfcEaJ^ 

Now, we introduce (on a possibly enlarged probability space where the conditional expec- 
tations Efc_i, Ejfc, Eaj. keep their definitions as the conditional expectations with respect to 
the same a-algebras as above) an x random symmetric matrix 

(56) A' = [a[j]i<ij<N 

such that: 

• the a'jj's such that i > k and j > k are i.i.d. copies of an (modulo the fact that A' 
is symmetric), independent of A, 

• for all other pairs (i, j), a'-j = aij, 
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then we have 

where is defined out of A' in the same way as is defined out of A in f l55|) (note that 
the k-th column is the same in A and in A'). It follows that 

NEk-i{YkYl) = EajEfc(/fcX/;;)]-EfcEaJfcxEfeE,J^ 

(57) = Efc[E,,.(/fc X /;;)] - EfcEaJfc X E^E.J^. 

We shall in the sequel prove that as N tends to infinity, regardless to the value of k, we 
have the almost sure convergences 

(58) E,Jfe L{z) E^Jl L{z') 

and that for any fixed u G (0, 1), as A^, k tend to infinity in such a way that k/N — > u, 
we have the almost sure convergence 

(59) ^.Mk^fk) ^ ^"(^,^0, 

where L{z, z') = \E'"(2, z'). The convergences of (l58l) and (l59i) are based on an abstract 
convergence result stated in next section, where we use the convergence of generalized 
moments of Proposition 12.11 They are stated in Lemmas 13.41 and 13.51 respectively. 

Note that by Lemma 16. 6^ 

(60) \fk\<M^z\-\ 

so once fl58l) and fl59|) proved, the convergences will also hold in L^, hence by continuity of 
the conditional expectations, we will have the convergence in 

Ek[E^Mkiu)y< fkiu))]-^kE.JkX^kE.J'k <iJ''iz,z')-L{z)Liz'). 

Thus by (1571) . we will have proved the convergence of (1541) . hence completing the proof 
of the theorem. 



3.3. An abstract convergence result. Remember that Ak = {dij) is the square matrix 
of size — 1 obtained by removing the /c-th row and the /c-th column of A, that Gk{z) = 
{z — Ak)~^ and that A'^ is a copy of Ak where the entries for i,j>k are independent 
of Ak and the other are those of A^. We denote G'i^{z') = {z' — A'f,)"^. 

Lemma 3.2. ^45 goes to infinity and tends to u in (0,1), the random probability 
measure on 

(61) ^z,? ■■= ]^3T hG,iz),„G',{z'),,}, z,z' eC\R. 

converges weakly almost surely to a deterministic probability measure ^, on C^. 
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Proof. By e.g. Theorem C.8 of [T], it is enough to prove that for any bounded and Lipschitz 



function /, i^'^'z'if) converges almost surely to i^zz'if)- Moreover, adapting the proof of 
Lemma [6 .4^ one can easily see that for any such /, 

converges almost surely to zero. The only modification of the proof is to complete the 
resolvent identity by noticing that 

{z-A,)-'-{z-A^)-' = {z-A^)-'{z{A,-A^)-{A,-A^r-2A^{A,-A^)){z-A,)-\ 

which gives that this matrix has rank bounded by 3 x rank(y4 — B). Hence it is enough to 
prove that the deterministic sequence ^[i^^'z' ( ' )] converges weakly. These measures are 
supported by the compact subset B{0, \Qz\~^Y, hence the sequence is tight and it suffices 
to prove that it admits at most one accumulation point. 

Let us show how the Levy case can be deduced from the case of Wigner matrices with 
exploding moments. For any B > 0, let us introduce the truncated matrices 

A^ := {aijl\a,j\<B)i,j=i,...,N A'^ := {a[ jliay<B)i,j=i,...,N- 

These matrices are Wigner matrices with exploding moments with rescaled moments that 
satisfy (jSj) (see [3T], Section L2.1] or [lOl Section 9]). Introduce the random probability 
measure i^^'^'^ , defined as i^^'j^, except than A and A' are replaced by A^ and A'^ , and 
similarly the matrices Gk.B{z) and G'f.Q{z'). By Lemma 2.4 of [TU], we know that for any 
e > 0, there exist B{e) > and, for any B > B{e), 6{e) > such that with probability at 
least 1 - e-''^^'^)^, 

rank(Afc - ) + Tank{A'^ - A'/) < e{N - 1). 

We can adapt arguments of [TU] (see also Equation (91) of [12] ) to get that with probability 
at least 1 — g-'^C^-^)^^ for any 1-Lipschitz function / on -8(0, IQ^^I"^)^, we have 

Hence it is enough to prove the weak convergence ofE[i/^^;^^( ■ )] when A is a Wigner matrix 
with exploding moments. Since these measures are uniformly compactly supported, it is 
sufficient to prove the convergence in moments. 

We consider a polynomial P = x^^xi'^^Xg^xi'™'^ and remark that 



E[<';(p)] = E 



1 ^-1 

^ ^ 5: (G, (z),,0 "HG. (.),,) (GU^O..) (GU^-')..) 



E 



— Tr |Gfc(z)°"i o Gfc(z)°'"i o G'fc(z')°"' ° G'^izT"^ 



_N- 1 

where o denotes the Hadamard (entry-wise) product of matrices and 

M°" :=Mo...oM. 
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We set ii = rii + mi and £2 = ^2 + m2. Let (YJ, l^')«=i,---/ij=i,---/2 be a family of random 
variables such that for any Pi,qj > 0, 

(62) Ejni^^^Hi;"' 



'Qj 

i=l j=l 



E 



N-1 



o . . . o Al!' O A'j^"' o . . . o A','* 



"=2 



Such a family exists by Proposition 16.81 By pT| Proposition 3.10], the couple of random 
matrices {A, A') satisfies the so-called convergence in distribution of traffics, which implies 
the convergence in moments of {Yi,Yj)ij. For reader's convenience, we give the limiting 
value of (162|) . even if we do not use it later. It is obtained by applying the rule of the 
so-called traffic freeness in 131). 



n 

i=i j=i TT&r(VT) 



(63) E[n^''n^"'i E at-mt- 



where 



1. we have considered T the graph whose edges are labelled by indeterminates a and 
a', obtained by 

• considering the disjoint union of the graphs with vertices 1, . . . ,pi and edges 
{1,2}, ...,{pi- l,pi},{pi, 1} labelled a, z = 1, . . . , £1, 

• considering the disjoint union of the graphs with vertices 1, . . . , and edges 
{1, 2}, . . . , {qj - 1, Qj}, {qj, 1} labelled a', j = 1, . . . , £2, 

• identifying the vertex 1 of each of these graphs (we get a connected graph, 
bouquet of cycles), 

2. we have denoted by Vt is the set of vertices of T, V{Vt) is the set of partitions of 
Vt and T'^ denotes the graph obtained by identifying the vertices of T that belong 
to a same block of vr, 

3. the quantity r°[T'^] is as in Proposition l2.lt 

4. we have set 



-"-/the edges hnking adjacent vertices'! W \^ 
1/ _T/ n/ 1 ^2 have the same label / 



V^=ViUV2 



where the sum is over all partitions of the set V.^ of vertices of T'^ and \Ei\ is the 
number of edges between adjacent vertices of V^, i = 1,2. 

Formula ( 163|) could also be derived by the same techniques than those developed in 
Section |2j The random variables Yi and Yj are distributed according to the limiting 
eigenvalues distribution of A. Moreover, following the proof of [381 Proposition 10], the 
exponential power series of this measure has a positive radius of convergence. So, by 
a generalization of [TTl Theorem 30.1] to the multi-dimensional case, we get that the 
distribution of {Yi, Yj) is characterized by its moments. Then, we get that {Yi, Y-) converges 



23 



weakly to a family of random variables {yi, y'j). We set fz{y) = {z — y) ^. We then obtain 
the convergence 



EK''(P)] = E 



i=i j=i 



E 



i=i j=i 



Hence E[i 



converges weakly. 



□ 



This convergence could also be proven without Proposition 16.81 but with appropriate 
bounds on the growth of moments. We have the following Corollary. 



Corollary 3.3. 



'!) For z E 

N-l 



, t SO that '^zt > 0, 



N 



^ ^{tGk{z)ii) pz{t) := J ^{tx)ul^,{dx, dx) a.s. 



(2) For z, z' e C\M, t, t' so that A = {f^z > 0, t'^z' > 0, |t| + \t'\ > 0} and k/N going 
to u, 

N-l 

N,k/ 

I,. I, I : = 

N 



converges almost surely towards 



i=l 



(64) 



^{tx + t'x')u^^,{dx,dx') . 



(3) For any {t,t') G M, the functions p^'^{t.,t') are analytic on A = {z,z' : f^z > 
0,t"^z' > 0} and uniformly bounded on compact subsets of A. They have a non 
positive real part. Their limits are also analytic on A and have a non positive real 
part. 

(4) For all z e C+ , 



1 r 

hm ^r{iz-Ar') = ^ / 



a.s. 



Proof. Let us first notice that p^(t) = p^f,(t,0) for all k so we only need to focus on 
P^At,t'). The point wise convergence of the function p^''^ is a direct consequence of the 
continuity of $ (recall we assumed that $ extends continuously to the real line), of the 
boundedness of tG k{z)ii+t' G'^{z') ^ (by t/Qz+f /Qz') and Lemma IX^ To show analyticity, 
note that for all j G {1, . . . ,N}, G{z)jj is an analytic function on C"*^, taking its values 
in C~ (and vice versa) almost surely by (l95l) . Hence, on A, tGk{z)ii + t'G'j^{z')ii is an 
analytic function with values in C~. Therefore, as $ is analytic on C~, p^'^ is an analytic 
function on A almost surely. Moreover, as $ extends continuously to the real line, it is 
uniformly bounded on compact subsets of C~ and hence p^'^ is uniformly bounded on 
compact subsets of A. This implies by Montel's theorem that the limit of p^'^ is also 
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analytic. Finally, p^'^ as non positive real part as the image of C by $. Indeed, as 
3?(— iA|aiip) < 0, we have 



(65) 



>N 



(A) I = |E(e-*^l""l')| < 1 m{\) = lim A^log |07v(A)| < 0. 



For the last point, first note that by Lemma 16.31 it is enough to prove the result for 
E[N-^TtG{z)] instead of N-^TtG{z). Second, by exchangeability, E[N-^TtG{z)] = 
E,[G{z)ii]. Remind that is the k-th column of A (or A') where the k-th entry has been 
removed. Using Schur complement formula and getting rid of the off diagonal terms (by 
arguments similar to those of Section [3l2|) . we have for any z G C"*^ 

E[G{z)n] = E 
= E 

Remember that we have for A G 
(66) 



^ - T.f=i \MjWGi{z)jj + En 



oil). 



1 

A 



sgn;^ oo 



and that 0Ar(A) = E[exp(— zA|aiip)] , which gives 



POO 

E[G{z)n] = I \ E[e*^(""^f="i'"^(^')'^^(")^-^)]dA + o(l) 



Af-l 



3*"^e[ J]0;v(AGi(z)j,)]dA + o(l). 



We used the exponential decay to switch the integral and the expectation. This also allows 
to truncate the integral: for any M > 0, 



N-\ 

e^^^E[n0^(AGi(^),,) 



dA 



M 



e^^^E[n0^(AGi(.),,) 



dA + e 



m,M,z,N-i 



where eM,z,N goes to zero as M tends to infinity, uniformly on and on the randomness. 
Remember that by assumption ([8]), we have 



(67) 



nUnW - 1) — y $(A), VA G c- 



where the convergence is uniform on compact subsets of C~. Hence, since for \'^z\ > S > 0, 
\t\ < M, tGk{z)jj belongs to the compact set {A G ; |A| < M/6}, we have 



N-l ^ N-l 

n MtGkizW = exp (- ^ HtG.iz),,)) + el% = exp(pf (t)) 
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where s^^l ^ converges almost surely to zero as goes to infinity. By Corollary 13 ■3[ 



't,z,N 

N-1 



Yl '^(tGk{z)jj)^ = exp (p,(t)) 



(2) 

where converges to zero almost surely. Hence, we deduce the almost sure convergence 
(68) e^*^ n MtG,{z),,) e^* 



7V-1 



As has non positive real part, we conclude by dominated convergence theorem and by 
getting rid of the truncation of the integral. □ 

3.4. Computation and convergence of Eaj./fe and Eafe/fc- 
Lemma 3.4. Almost surely, we have the convergence 

(69) J, Liz) := - / ^.e^^^+^^-^dt . 

Af->oo Jo t 

where sgn^ := sgn(Q'2;) and Pz{t) is defined in Corollary \3.3[ 
Proof. Remember that for any z G C\M, 

^ z -T.j^k{jyGk{z)jj ■ Xn{z) 

Hence, by ( 166|) . since by (p5|) the sign of the imaginary part of Xn{z) is sgn^, the random 
variable fk can be written 

/•sgn^ oo 

(70) /, = -W 9,A;v(^)e'*'"(^)dt. 



^0 

= -z / 9,A7v(z)e'*^^(")dt - z / dANiz)e''^''^'Mt 

J sgn^ m ^0 
• fk,m ~l~ Vmi^z) J 

where m > and r/^(z) = - e^'^g"-'^^^^^)). We next show that for all e > there 

exists mo small enough so that for all m < mo, all large enough 

(71) ^.A\VmM]<e 

By (197|) and since the sign of the imaginary part of \n{z) is sgn^, one has |77m(2^)| < 4|32|~^. 
More precisely, for any > 0, we find 

4 
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In the exploding moment case, we note that 



K 



whereas in the Levy matrices case, since Yl^kUY converges weakly towards a a/2 stable 
law, we know that for large enough 

C 



Choosing K = 1/y/rn or K = m 2+a we conclude that there exists a 7 > and a finite 
constant C which only depends on '^z > so that for N large enough 

(72) E^,[\vM\] < Cm'' . 

In the following we shall therefore neglect the term rim{z). By f l^^ and f l9B]) we have the 
following bound: for any t 7^ and any 2; 6 C \ M such that \'^z\ > 5 > 0, 

(73) I afc(jr(G,(2)2),.,.e-*S^'=(^-)'^'=(^)-- 
so that for M large enough and \'^z\ > 5 > 0, 



< - sup xe 

x>0 



-tx 



ts 



—I 



sgn^ M 

where em,M,z,N is arbitrary small as M is large, uniformly on N and on the randomness. 
Moreover, by f l75]) one has 



(74) 



Recall that </)Ar(A) = E[exp(— zA|aiip)] , so we have 



(75) 



Ea, [e**^-(^T = Ea, 



Atz 



N-l 



Remind that in fl68il we have shown in the proof of the last item of Corollary 13.31 the 
almost sure convergence 

Af-l 

(76) (^'^''Y[^.r(+nJA..\ ^ ^itz+p,{t) 



e^*^ n MtG,{z),i) 



As in the proof of Corollary 13. 3 [ since the left hand side is analytic and uniformly bounded, 
we deduce by Montel's theorem that its convergence entails the convergence of its deriva- 
tives. We then get by fl75|) . for all t, z so that tf^z > 0, the almost sure convergence 



dJ E 



N-*oo 
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We then obtain by dominated convergence (remember the integrant is uniformly bounded 
by dZH) and ((73])), 

f sgn, M 

'm,M,z,Nj 

T 

I sgn^ m 

where em,M,z,N converges to zero almost surely as goes to infinity. 



/•sgn. At -| 

^.Jk,m = - ~ -S^e^^^+^^Wdt + e, 



By f lTIjl and estimate fl73|) . we have the estimate 



so we can let M going to infinity to obtain the almost sure convergence 

^.Jm,k ^ - / -9,e^*^+''^Wdt. 

In the Levy case, one has Pz{t) = Pz{^), and in the exploding moments case, | jS^e**^"^''^''*-' | < 
1 + so that the integral converges at zero and we obtain 

Af^oo Jo t 

as expected. □ 

Of course, this convergence is uniform in k since the law of Gk{z) does not depend on k 
and an analogous formula is true for Eafe/fc, replacing 2; by 2;'. 

3.5. Computation of Ea^,(/fc x fj!). 

Lemma 3.5. Almost surely, we have the convergence 

/•sgn^ 00 /"Sgn^/ 00 -1 

(77) EaJ./^' vl;«(,,,'):= / / i;a,Ve'"+'*'^'+<-(*'*'^dtdt', 

N^oo Jq Jq tt' 

where p"^, is defined in Corollary \3.3[ 

Proof. We shall start again by using Formula f lTOj) for f^, and its analogue for 

. ^ l + E,a.(jr(GU^O%- ^ g.A-^(zO 
^-E,a.(jrGU^').. ■ A'^(2') ' 
where ^^(2') is defined as Gk{z), replacing z by z' and the matrix A by the matrix A' 



defined by (l56l) . which gives 

/■sgn^, 00 

(78) = -I 9,A'^(2')e^*'^-(^')dt' 



— ? 





sgn!, 00 /"Sgn . m 

dzX'^{zy'^'^^'"^dt' -I / ^ dz\'^{z')e'''^'^^''^dt' 

sgn', m ^0 



28 



FLORENT BENAYCH-GEORGES, ALICE GUIONNET, CAMILLE MALE 



The upper bound (1731) allows us to bound the first term uniformly by logm^^ and to 
truncate the integrals for sgn^ t, sgn^, t' < M. Therefore, up to a small error e uniform for 
M large, m small and provided I'^zl, \ Qz'\ > 5 > 0, we have 

J sgn^ m J sgn^ m 

As in the previous case, the upper bound fl73p allows us to write 

1 



(80) E; 



tt 



Remember that 0Ar(A) = E[exp(— zA|aiip)] , so we have 



(81) 



E, 



JtXNiz)+it'\'{z') 



e^'^+^''^'l[<j>4tG,{z),,+t'Gi{z%). 



By assumption ([HD and Corollary 13.31 we have the following almost sure convergence as 
goes to infinity and 4 goes to m in (0, 1) 



N~l 



^tz+u'z' -Q ^^(^tG,{z),,+t'G[{z%) e 



itz+it'z'+p^ _,(t,t') 



Almost surely, for any t,t', the map {z,z') t-)- e**^"*"**'^' (^G^fc(^)ij + ^'^U^Oii) is 
analytic on C^^"* x <C^s^^' and bounded by one. Hence, with the same arguments as in the 
previous section, we get almost surely and for any t, t' so that t/Qz and t' f^z' are positive, 
the uniform convergence for the second derivatives on compact subsets. The truncations 
of the integrals can be suppressed as in the previous section, we obtain the almost sure 
convergence 

fsgri- oo psgn^i oo 



N-^oo 



dtdt'. 



□ 



Hence we have proved the convergences 058!) and flSUj) . This completes the proof of 
Theorem 11.31 



4. Proof of Corollary 11.4! 

By linearity and the fact that A is stable under complex conjugation, it is enough to 
consider a single function f & A. By Cauchy formula and Hypothesis (ITOi) on /, we have 
for all X G M, 



29 

where F = M ± ie. Let = /rn{|Kz|>L} fS'^'^ ^^"^ observe that for L big enough, 

the polynomial decay of / at infinity gives 



Therefore, Lemma [6.31 shows that ZNifi) goes to zero in probability as TV goes to infinity 
and then L goes to infinity: for any rj > 0, 

lim limsupP(|Z^(/L)| > rj) = 0. 

L-s>oo N^oo 

Fix L,5 > and consider Xk = —L + k6 for G [0, i?], with i? = 12LS~^\. Let 

217, J. g^.. z-x 2z7r^f^^^A^ 

Note that for each L and e > there exists a finite constant e) so that 

\\h-Il\\TV<K{L,e)5. 

Hence, Lemma 16^31 implies that Zn^/l) — ZN{fl) goes to zero in probability as goes 
to infinity and S goes to zero for any given L, e > 0. Finally, by Theorem ll.3[ ZN{f[) 
converges in law towards a centered Gaussian variable with covariance 

f(y + '>'Se)dy / f{y + is'e)dyC{xk + ie, Xk' + ie) 

The above covariance converges as 6 goes to zero by smoothness of the functions C towards 

-— E] / / dxdx f{x + ise)f{x' + is'e)C{x + zse, x' + is'e) . 
^^^,^^J\x\<lJ\x'\<l 

Therefore, writing 

= + [ZmUl) - Zm{R)] + Zn{R) 

and fixing L large enough and then 5 small enough so that the first two terms are small 
with high probability, we conclude that ZnU) is close to ZnWl) and therefore converges in 
law towards the announced centered Gaussian variable. The generalization to the closure 
of A for the total variation norm is a direct consequence of Lemma 16.31 



5. Fixed point characterizations 

In this section, we provide characterizations of the functions pz and ^, involved in the 
covariance of the limiting process of Theorem 11.31 as fixed points of certain functions. Levy 
matrices and Wigner matrices with exploding moments are considered in a unified way. 
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5.1. Fixed point characterization of p^( ): proof of Theorem 11.51 We now prove 
the fixed point equation for the non random function involved in the Lemma 13^ given for 
2; e C+ and A > by, 



p,(A)= lim p^^(A) 



N 



hm 



1 ^ 

-5^<|.(AG.(2),,), 



where we have proved that this convergence holds almost surely in Corollary 13. 3[ Note 
however that under the assumptions of Theorem 1 1.5 1 the arguments below provide another 
proof of this convergence. 

We denote in short A for A^, G for Gk and a for a^ in the following, and we do not detail 
the steps of the proof, which are very similar to those in [TU] and Corollary 13.31 but outline 
them. Since we have already seen that p is analytic in Corollary 13.31 we need only to prove 
the fixed point equation. Let Ai be the N — 2 x N — 2 principal submatrix of A obtained 
by removing the first row and the first column of A, and let Gi{z) := {z — Ai)"^. Let ai be 
the first column of A where the first entry has been removed. Using first the concentration 
lemma [631 then exchangeability of the G{z)jj^s, then Schur complement formula (see [H 
Lem. 2.4.6]), and then the fact that we can get rid of the off diagonal terms by the same 
argument as in the proof of Lemma [XT] since $ is continuous on C^, we have for all z G 

E[pf(A)] = E[<|.(AG(z)n)] 



E[$( 



A 



an 



E 



a*Gi(2;)ai 
A 



)] 



n 



0(1). 



Then we use Hypothesis made at (1181) to get for '^z > 



E[pf(A)] 



^f{^-E,ai(j)2Gi{^),,) 



]dy + o{l) 



V 

A^E 



N-l 

n 



'^N{jGi{z)jj] 



dy + o{l) 



Using the definition of $ and the fact that we assumed that it is bounded on every compact 
subset (since has non positive real part we can cut the integral to keep y bounded up 
to a small error, as in the previous sections), we have 



E[pf(A)] 



g{y)e'-^E 
9{y)e 



dy 



dy + o{l) 

0(1) 



dy + 0(1) 
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Now, since 3f?(p^~^) < and g satisfies (fT9|) . we may use Corollary 13.31 to take the large 
N limit and we get ( 120|) by dominated convergence theorem (by Lemma \Q.2\ p^~^ can be 
replaced by in the above formula). 

Let us now prove that there is a uniquene solution to this equation which is analytic in 
C+ and with non positive real part. Suppose that there are two such solutions p^(A), p^(A). 
For z fixed, let us define A(A) := \pz{\) — pz{X)\. Then for all A, we have, by the hypothesis 



made on g in ( 1T8|) . 

POO 

A(A) < A / \gi\y)\e-y^^Aiy)dy 
Jo 



It 



< K\^X^'+^ j y^'e^y'^'A{y)dy + X^+^ j y''e-y'^'A{y)dy 

(poo POO \ 

A^+^y^ y^e-y^^A{y)dy + X^~^' y''e-y^^A{y)dy] . 

POO POO 

follows that Ji := / A^e-^^^A(A)dA and h := / A^e"^^^A(A)dA satisfy 
Jo Jo 

(POO POO \ 

Ji J A^^+^e-^^MA + I2J A^+^+^e-^^MA J , 

(POO POO \ 

h J A^+^+^e-^^MA + I2J A^'^+^e-^^MA J . 

For large enough, the integrals above are strictly less that 2^), so for large enough, 
Ii = I2 = 0. It follows that for any fixed A, p^(A) and Pz{X) are analytic functions of z 
which coincide for large enough, hence are equal. 

5.2. Fixed point characterization of Pz,z'{-, ■)'• proof of Theorem II. 6L We now find 
a fixed point system of equations for the non random function of Corollary 13.31 For 
X/Qz + X'/Qz' > 0, we set 

(82) Pz,t'\K^') = ^Y.^{XG,{z),, + X'Gi{zX, 



j=k 



where we recall that Gk and G"^. are as in f[T^ . To simplify the notations below, as in 
the previous section, we denote (G, G') instead of {Gk, G'l^), even though their distribution 
depends on k. 



In the sequel we fix, as in Section 13. 3[ a number u G (0, 1) and will give limits in 
the regime where — ?■ 00, A; — 00 and k/N — y u. We shall then prove that, under 
the hypotheses of Theorem 11.61 that we assume throughout this section, {pzV^ ^ PzV"^) 
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converges almost surely and that its limit satisfies a fixed point system of equations which 
has a unique analytic solution with non positive real part. The convergence could be 
shown with minor modifications of Lemma [3l2l but we do not need this since we work with 
stronger assumptions. Using the concentration lemma 16.41 (note that $ is not Lipschitz 
but can be approximated by Lipschitz functions uniformly on compacts), it is sufficient 
to prove the fixed point equation for the expectation of these parameters. Moreover, by 
exchangeability of the k first entries and N — k last entries 

^Pz}\t,s)] = E[<l>{tG{z)^M + sG'{z')NN)]+o{l). 



These functions are analytic in A = {z,z' : t/Qz + sf^z' > 0} and uniformly bounded 
continuous in {t,s) (by uniform continuity of $ and boundedness of G{z)), and hence 
tight by Arzela-Ascoli on compacts of A. We let p^'li,i = 1,2 be a limit point. Since 
K[p^p\t, s)'\ is uniformly bounded on compacts of A, the limit points p"'*,z = 1,2 are 
analytic by Montel's theorem. We assume for simplicity that both z and z' have pos- 
itive imaginary parts (in the general case, one only has to replace /t=o° below by 
^^sgn(s.)oo^^sgn(s.')oo^_ Under Hypoth esis made at (1241) . we can write by Schur complement 
formula and getting rid of the off diagonal terms (note that all integrals are finite as they 
contain exponentially decreasing terms) 



oo poo 



JO 



e*iE[e"''i^^"i(^)'^i(")«]d/i(tO + J e'iE[e-'^^^"^(^)'^'=("')«]d/i(i;) + o(l) 



oo roc 



JO 



POO 

+ / e'?e''^(?M/i(i;) + / e'^e^^'(?M/i(^;) + o(l) 



Jo 



and 



E[p^;^2(t,s)] = / / e*f^+*T-'E[e-*i^^^'-W'^^(^)"-^^^^^^W'^'^(^)«]dr(i;,i;') 
Jo Jo 

+ J e*TE[e-'?^^"^(^)'^'^(")«]d/i(i;) + J e*?E[e-'^^^"'^(^)'^*("')«]d/^(w) + o(l) 



oo roo 



JO 

oo 



+ / e'^-e'^^^-^dfi{v)+ / e'?e^^'(?Mp(?;) + o(l) 
Jo Jo 
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where we used that ai(^) = Si[{£) for all £ whereas aAr(£) = a'^(£) for i < k and are 
independent otherwise, that converges towards pz, and that 



k-2 

k-2 



N,k,2 



j=k 



by Lemma [6.21 (bv continuity of $, it can be approximated by Lipschitz functions). Hence 
we find that the limit points p"']/,^"^/ of p^'^i'^ , P^t''"^ satisfy ( 126|) . Moreover, 



pz:z' = upzjz' + (1 - + 

gives 

n M.l I /-I \ u,2 

Pz,z' = UPz]z' + (1 - u)Pz\z' ■ 

Uniqueness under assumption (1): Let p and p be solutions of Equation f l26p with non 
positive real parts (note here that Pziu) is given and /i is so that the above integrals are 
finite; hence the last two terms in both equations play the role of a finite given function). 



fOO POO 



< 2/ / e----^--^'^'-^^*l^A(^,^)d.d.' 



Jo 

oo POO 



POO POO 

< 2ts / / e-^^''-^^'^'i^((^;t)^l,te]o,i] + (^^t)"l.te]i,oo[) 
Jo Jo 

x[{v'syi.y>se]o,i] + (f's)''lj,'s6]i,oo[)A(t',t;')dtdt/ 
where for (a, a') G {7, k}^, we have set 

POO POO 

lAia,a'):= / e-'''''-^''''' Aiv^v'^'^v'"' dvdv' . 
Jo Jo 

We put I {a, a') = Ii{a,a') where 1 denote the constant function equal to one. We get 
after integrating both sides 

lA{a,a') < i^(/(a + 7 + l,"' + y + 1)^a(7,7) + ^(" + 7 + + + 1)^a(7,^) + 
I{a + K + 1, a' + 7' + 1)/a(k, -f) + I{a + K + l,a' + k + 1)Ia{i^, «:)) • 

We consider Qz, ^z' large enough so that a'+f3'+l) < for any (a, a'), /?') G 

{7, k}^ to conclude that A(t, s) vanishes then and therefore that p"^ = p"'^, for 0^2; and 
Q'z' big enough, s = 1 or 2. By analyticity, we conclude that the system of equations ( 126|) 
has a unique analytic solution with non positive real part. The functions p^''^'^ and p^^^'"^ 
are tight and their limits points are characterized by fixed point equations, so they actually 
converge. 
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Uniqueness under assumption (2): We next consider the case where = cx". 

The hmit points of {p^'z''^, P^f/'^) satisfy the system fl26|) with r given by (128|) . that is 



dr(f , v') 
dvdv' 



+00 







To simplify the notations we assume hereafter ^z' non negative. Notice first that 
if ((71, (72) is a hmit point of {p^p^ , P^'z''^), Qi are function from iW^Y into Ca/2 '■ = 
|_^g«f ^ < a/2}, i = 1,2, which Qi are homogeneous of degree ^, i.e. for any t,s > 

t s 



gi{t,s) = {t' + s')^gi 



We show that for any 13 in (|, 1), the system (126|) has a unique solution on the set of 
pair of homogeneous maps (M"*")^— )'£a/2 of degree | that satisfy the /9-Holder properties, 
i.e. have finite || ■ norm given by 



^5) 



1^1 



max \g{u, v)\ + 



max 



iu,v)GSl 

with Si = {s,t> 0,5^ + ^2 = 1}. 





-9 


yy)\ 




-{v 


- u')2|/3/2 



Proving that the limit points (p"^/, p"^/) of (p^^^^, P^fr'"^) are /3-Holder maps for a /3 > | 
allows to conclude the proof of Theorem 11.61 This is the content of the following lemma. 

Lemma 5.1. For any z, z' G C+, m G [0, 1], i G {1, 2} and f3 G {a/2, (3a/4) A 1), 

II u,i II . 



Proof. First, since |G'(2;)«| < '^z ^, we have 



(86) 



max max 
i=i,2 (s,t)eS'i 



<l^+ ' 



'^z '^z' 



We next show that for any matrix model corresponding to $(a;) = — (t(zx) 2 we have for 
2k G (0,a/2) 



^7) 



Indeed, we can write 



limsupE ( 



N>1 



2k 



< 00 . 







= cE 













-yT,\o-ii? 



,1+2k 



-dy 



-dy 



where we have used Fubini for non negative functions. But the above integral is well 
converging at infinity and we know that (p^ converges uniformly on [0, M] for all M finite; 
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hence there exists a finite constant C so that for large enough 



sup |y|-tiV|0^(-2y)-l| <C 

j/6[0,M] 



which yields 



^^^^ dy ) < oo for 2k < a/2 



'0 y 

We next show that this estimate implies the /3- Holder property. Indeed, for any /3 G [f , 1], 
there exists a constant c = c(a, /3) such that for any a;, y in C~, 

(88) ~V^\ <c\x- yf (|x| A \y\)^~'^ ■ 

Applying this with x = tG{z)jj + sG'{z')jj and y = t'G{z)jj + s'G'{z)jj gives 



< cK 



N 



+ 



3^2 (CJ^/)2 



U\t-t'\' 



\S — S 



f) 



with 



Km := E 



{\tGiz)^^ + sG'iz')-n\ A |t'G(z)n + s'G'iz')^,\y 



where k = /3 — |>0. It is enough to prove that Kn is uniformly bounded as p"'*, is a 
limit point of E[p^i,^'*]. Note that we may assume that |s — s'| < 1/6 and |t — t'| < 1/6 
since otherwise the bound is already obtained by fISB]) . But then this implies that either 
t,t' > 1/4 OT s, s' > 1/4. Let us assume t,t' > 1/4. Then, we have 

\t'G{z)u + s'G'{z')u\ A \tG{z)u + sG'{z')ii\ >t' A t^G{z)u > ^\QG{z)u\ ■ 
Using Schur formula, we find that 

- a*Gi(2)ai))2 + {^{z - an - a*G'l(^)al))2^ ^ 



EpGiz)ur] = E 



< 



2K 



E 



'^{z — a^Gi(2;)ai^ 



so that we deduce that for all k > 



E[|^>G'(^)nr''] < 



E 



Equation (187|) completes the proof by taking 2^ = 2(/3 — f ) < f • 



□ 
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We now prove the uniqueness of the solution of (126|) on the set of functions described 
above. After some change of variables, the system is equivalent to the following: 



fOO PCO POO 



Jo Jo Jo 



foo poo roo 



g,{t,s) = / / y;«/2-l(^')«/2-l^-a/2-le.l^.^ 

Jo Jo Jo 



ivT.Z+ugi{W+vT)+{l-u)p^_^,{W+vT) _ e«9i(W")+(l-")Pz,.'(W^) ) dt^'df 



where we have denoted in short W = {w,w'), Z = {z,z'), T = {s,t), Pz,z'(W) = Pziw) + 
Pz'iw') and W.Z,T.Z stand for the scalar products. 

After the change of variables w = r cos{6),w' = rsin(6'), v = rv', we can rewrite this 
system of equations as 

91 = K,z'iugi + (1 - u)g2) 92 = Fl^,{ugi + (1 - u)pz,z') 

with if T'^ = [t, s) when T = (s, t), 

F^MiT^) = [ r ^-/2-le-,.Z ^^rvT.Z+r'^/^9iee+vT) _ ^r'^/^9{ee)\ ^^^^^ 
Jv,e Jr=0 ^ ^ 

where we have denoted eg = (cos(6'), sin(6')) and 

dp{v,9) = C„le6[o,|]rf^(sin2^)t-^l^e[o,oo)^^"^"^ 

for a constant Ca- 

The desired uniqueness follows from the following lemma. 

Lemma 5.2. Let /3 in 1). For any M and for '^z, '^z' large enough, the map F^^, is 
a contraction mapping on the set Cm,I3 of homogenous maps g : (M"*")^— 7'£q/2 of degree f 
with (3 norm bounded by M . 

Proof. To study the Lipschitz property of F"^, as a function of g in Cm,i3 for the norm /3, 
we first set 

F:A9){T\f^) = FlAg){T^) - FlAg){f^) 

= j dp{9, V) j drr^/'^'^e"'""-^ ^^irvT.Z+r°'/^gieg+vT) _ ^irvf.Z+r°'/^g{eg+vf)^ 

We next bound, for given gi, g2 in Cm,/?, Ti in S]_ and with T2 either in S\, or T2 = (which 
allows to treat in one time two parts of || ■ ||^) 

Af" = KAgim,T2) - F:jg,)iT,,T2)\ . 
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For that task, we shall use some technical estimates, with a constant c that may change 
from line to line. Remind first the two following bounds from [121 Lem. 5.7] valid for 
7 > 0: there exists a constant c = c(a,7) > such that for all h,k & C"^ and for all 
xi,x2,yi,y2 e £s. 



(89) J'y,h{xi,X2,yi,y2) 







< c (|/ir^~^|xi -X2- yi + y2\ + |/ir^"°(|xi - X2I + \yi - y2\){\xi - yi\ + \x2 - ^2!)) 
Moreover, for all h,k & C"*", x, ?/ G £q , for < k < 1, we have 



K^,h,k{x,y) := 

(90) < c{\h\A\k\y"'~'i-''\h-k\''\x-y\, 

Moreover, notice that for {s,t) in S], one has max(s, t) > 1/72 and max(cos6', sin6') > 
1 / -\/2 and thus 

ST.Z >QzA Qz'/V2 =: A,,,/ '^{ee-Z) > A,^,,, 
so that \ie0.Z + zi;T.Z| > A^^^/(l + Ij^g^i u). At last, with Oj = + t;Tj for z = 1,2, 
straightforward uses of the /3-norm and the inequalities |aj| < (1 + f), ||ai| — |a2|| < 
^^iTi-Tal, |a,| > ^(^^ V 1) if T, e St, |^ - ^| < {vV iy^\T, - T2HI + v), and ([88]) 
gives the estimates 

2 

|5^(-l)^'+^-^7.(a,)| < c||^i-^2||,|Ti-r2|^/,(^), 

2 2 

( J2 - ( Yl - ^*(«2)|) < cM\\g2 - gih\T^ - T2f{l + v)'^fp{v), 

1=1 1=1 

where fpiv) = v^(^(l + v)^~^'^{v V + 1^ Using this series of estimates, we find that. 



Apu < / K 



^F" ^ / J:^^,eg.z+vT.z,eg.z+vf.zi9{ee + vT), 



vf))dfi{v,9) 



J<^^eg.z+vT.z{g{e0 + vT),g{ee + vT),g{ee + vf),g{eg + vT))dfi{v,9) 
< c\\g-~gy(^A:^^,lf^,l h>^dfi{9,v) 

+\f- TfA;;:r^{\z\ + l^'l)'^ |(1 + .;l|^|^,)-^-"(l + if=o,.<iyMo, v) 
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While using ( 190|) . we chose = ajl and k = 13 when Tg 5*1 or T = 0,f<l,/i; = when 
T = and v > 1. As the integrals are finite we obtain the desired bound for g,g & Cm,p 

WKAai) - FlA92)h < C{z,z',M)\\g-~g\\p, 
with C{z, z' , M) <1 ii^z /\ "^z' is large enough. □ 

Taking two solutions of the system (l26i) in Cm,!^, we deduce that they are equal when 
'^z A '^z' is large enough, and thus everywhere by analyticity. 



6. Appendix 



6.1. Concentration of random matrices with independent rows and linear alge- 
bra lemmas. This section is mostly a reminder of results from |12j and [T3] . 

The total variation norm of / : M — ?■ M is 

I TV := sup^ - f{xk)l 



where the supremum runs over all sequences {xk)k£Z such that Xk+i > Xk for any k E Z. 
If / = l(_oo,s] for some real s then ||/||tv = 1? while if / has a derivative in L^(]R), we get 



TV 



\f{t)\dt. 



The next lemma is an easy consequence of Cauchy-Weyl interlacing Theorem. It is an 
ingredient of the proof of Lemma 16. 3[ 

Lemma 6.1 (Interlacing of eigenvalues). Let A be an n x n hermitian matrix and B a 
principal minor of A. Then for any / : M — )■ C such that H/Htv < 1 one? lim|2;|_!.oo f{x) = 0, 



n-l 



i=l 



i=l 



< 1. 



Lemma 6.2. Let Ai, A2 be n x n random Hermitian matrices and Ai, A2 ben — Ixn — 1 
matrices obtained from Ai and A2 respectively by removing the i-th row and column, for 
some £ G {1, . . . , n}. Let z, z' G C, t, t' G M so that ^zt > and '^z't' > and set G = 
{z - Ai)-\ G' = {z' - A2)-^ and G 
f on B^^^>^t,t' 
have 



z — Ai) ^,G' = {z' — A2) ^. Then, for any function 



{g eC^ ; \g\ < G{z, z' ,t,t')} with G{z,z',t,t') = tC^z)'^ + t'{Qz' ) % we 



(91) 



k=l 



n 



k=l 



< 



G{z,z',t,t') 



n 



Lip 



n 



where WfWup ■= sup^^^ ^^^^z^p^ and 
the elements of Bz,z',t,t' ■ 



:= sup^ 1/(3;) I, both supremums running over 
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Proof. The proof is similar to [T2| (91)]. We denote by Ai,A2 the n x n matrices whose 
entries are the same as Ai, A2 except for the ith rows and column which have zero entries. 
We denote G, G' the corresponding Stieltjes transform. Then, G, G' equal G, G' except 
at the ith column and row (where it is equal to ;z~^li=j=fc). Therefore, noting M = 
tG{z) + t'G'{z') and M = tG{z) + t'G'{z), we conclude that 



- J2 fiMkk) f(^^>'^ 



k=l 



k=l 



< 



n 



Moreover, let M = tG{z)+t'G'{z') and note thatAi — and A2 — A2 have rank one so that 
M—M has rank one. On the other hand it is bounded uniformly by C{z, z',t,t') = ^ + ^- 
Hence, we can write M — M = cuu* with a unit vector u and c bounded by G{z, z', t, t'). 
Therefore, since / is Lipschitz, 



1 " 



k=l 



^ n—1 
k=l 



< 



< 



1 " 

Lip~ / > l^kk — Mkk\ 
k=l 

1 " 

Lip- ^ C{Z, z', t, t') (Cfc, uf 



k=l 

C{z, z',t,t' 



Lip 



n 



Equation (19T|) follows 



□ 



Lemma 6.3 (Concentration for spectral measures [H]). Let A be an n x n random Her- 
mitian matrix. Let us assume that the vectors (Aj)i<j<„, where Ai := {A.ij)i<j<i G C*, are 
independent. Then for any / : R — )• C such that ||/||tv < 1 and E| f f d^A\ < 00, and 
every t > 0, 



P 



/ diiA f >t] < exp - 



nt' 



Lemma 6.4 (Concentration for the diagonal of the resolvent), a) Let A be an n x n 
random Hermitian matrix and consider its resolvent matrix G{z) = {A 



, z e 



Let us assume that the vectors (Aj)i<j<„, where Ai := (Aj)i<j<« ^ C*, are independent. 
Then for any / : C~ — >■ M such that ||/||Lip < 1, and every t >0, 



P 



-J2f(G{z),,)-E-J2fiGi^)kk) 



k=\ 



k=\ 



>t \ < 2 exp 



n^(z\H^ 



b) Let A' be annxn self-adjoint matrices so that A'^^ = Aij,iAj < k and {A'j^j)j>k+i^i>k+i 
is independent from (Aij)j>A;+i,j>fc+i but with the same distribution. Let G{z) = {z — A)~^ 
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and G'{z) = {z — A')~^ and set for a Lipschitz function f on C , 

£=1 

1 ^ 

e=k+i 

Then, for Xf^z > 0, \' f^z' > 0, we have for all 6 > , s ^ {0, 1}, 

<S^((fc-l)^~°+(iV-fc-l)") 

(92) p(|p^J'^+^(A,AO[/]-E[p^;'^+^(A,A')[/]| > 5) < 2e «ii/iiLp^(^-^^'.*')^ 
with 

CiX,X',z,z') = ^ + ^^. 

ksz ^z' 

Proof. The first point is proved as in [121 Lemma C.3]. We outline the proof of the second 
point which is very similar to [T^ Lemma C.3]. We concentrate on p^'y^i the other case 
being similar. By Azuma-Hoefding inequality, it is sufficient to show that 

X, :=E[p:;J'^(A,A')[/]|J,] -E[p:^;y(A,A')[/]|Jp_i] 

is uniformly bounded by ||/||LipC(A, A', t, t')/c~^. Here J-p is the a-algebra generated with 
respect to the p first column (and row) vectors. Note that Xp can be written as the 
conditional expectation of the difference of the parameter / evaluated at two sets A, A 
and A^A' which differ only at the p-th vector column (and row). Hence, we may follow 
the proof of Lemma 16.21 to conclude that 

LiPV^I/.^ I _ |c|ll/llLip ^ „ ^ ^2^ |c|||/l|Lip 



□ 



Let H = [hij] he an N X N Hermitian matrix and z G C\M. Define G := {z — H) 

Lemma 6.5 (Difference of traces of a matrix and its major submatrices) . Let Hk he the 

submatrix of H obtained by removing its k-th row and its k-th column and set Gk '■= 
{z — Hk)~^ ■ Let also be the k-th column of H where the k-th entry has been removed. 
Then 

(93) Tr(G') - Tr(Gfc) = + ^l^l^k 

z — — alGk^k 

Moreover, 

(94) I Tr(G) - Tr(Gfc)| < 7r|^5z|-\ 

Proof For ([93]), see [H Th. A.5]. For ([94]), see Lemma lOl □ 
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Lemma 6.6. With the notation introduced above the previous lemma, for each 1 < j < N , 
(95) 



X QGjj < 0, 



(96) 

and for any a = (ai, 
(97) 



C>-y|-l 



1 + |a.f (G^- 



< 2 



Proof. Set z = x + x, 1/ e M. Let Ai, . . . , Aat be the eigenvalues of H, associated with the 
orthonormahzed collection of eigenvectors ui, . . . , uat. Let also ej devote the jth vector of 
the canonical basis. Then (195|) and ( 196|) follow directly from the following: 



N 



|(e„Ufc)|2 A |(ej,Ufc)|2 



(98) \{G' 



_ \, ^2 I - _ 



N 



l(ei,Ufc)|2?/ 



(^-A,)2'-^(x-A,)2 + y2 n ^(^_A^)2 + ^2 

Let us now prove flWl) . By flUSl) . we know that 53z and —'^Gjj have the same sign, so 



(99) 



1 



< 



1 



Hence it remains only to prove ( 197|) . This is a direct consequence of (195|) and ( 196|) which 
imply the second and last inequality 



< 



Zlj |aiPl(G'^)iil 



< 



E,|a,f|(G%.| _ 1 



< 



53(^-E,|a,fG,,(^))| - E,|a,f|53G,,| " ' 



□ 



6.2. CLT for martingales. Let (J^t)fc>o be a filtration such that J-'o = {0, f^} and let 

(Mfc)fc>o be a square-integrable complex-valued martingale starting at zero with respect to 
this filtration. For /c > 1, we define the random variables 



and we also define 



fe-i 



k>l 



Let now everything depend on a parameter A^, so that J-^ = J^k{^)-!^k = Yk{N),v = 
v{N),T = T{N),L{e) = L{e,N),... 

Then we have the following theorem. It is proved in the real case at [HI Th. 35.12]. 
The complex case can be deduced noticing that for z E C, ^{z)'^, '^{zY and 3fJ(z)Sj(2;) are 



linear combinations of z"^, z"^, Izl"^. 
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Theorem 6.7. Suppose that for some constants f > 0, r G C, we have the convergence in 
probability 



and that for each e > 0, 



v{N) — > V t{N) — > T 

N^oo Af— s>oo 



L(e,N) — > 0. 



Then we have the convergence in distribution 



k>l 

\2\ _ „, TT7/"72\ 



where Z is a centered complex Gaussian variable such that E(|Z| ) = v and E(Z ) = r. 
6.3. On the Hadamard product of Hermitian matrices. 

Proposition 6.8. Let Ai,...,Ap be N by N Hermitian random matrices whose entries 
have all their moments. Then, there exists a family of random variables ( ) whose 

joint distribution is given by: 



E[a1' . . . a";"] = E 
where o denotes the Hadamard (entry-wise) product 



Vni , . . . , np > 



Proof. Step 1. We first assume that the matrices are deterministic and have distinct eigen- 
values. By the spectral decomposition, for j = we have Aj = 'Yl!i=i^j,i'^3,i^*j,i 
where Aj = {\j,i)i=i,...,N is the family of eigenvalues of Aj arranged in increasing order, and 
— {.'^j,i)i=i,...,N is the family of associated eigenvectors. For any ni, . . . ,np > 0, one has 

1 1 N 

-Tr[Aro---oA;^] = -$^(Ar)(A:,fc)...(A;^)(fc,A;) 

k=l 

N / N X /TV 



1 / \ / 

k=l ^n=l ^ ^ip=l 

ii,...,jp=l ^ k=l ^ 



For j = 1, . . . ,p, we set dyUA^ = Xlili ^A^i the empirical eigenvalues distributions of Ay 
We denote Fa^. it) = /iA^ ((~oo, t]) the cumulative function of d/iAj - Since the eigenvalues of 
the matrices are distinct, one has F\_.{Xj^i) = for any i = 1, . . . , N, j = 1, . . . ,p. Hence, 
we have 

i- Tr[Ar o ■ ■ ■ o A;^] = [ XT... X;^fN ((A„ A„ f/,),=i,...,p)d/iA, (Ai) . . . d/iA,(A,), 
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where /^((A„A„f/,),=i„„,,) = (^N^~'EtiUU\ 

%(A'i^Aj (Aj))(^)r^ • Hence, a family of 

random variables (ai,...,ap) as in the proposition exists and its joint distribution has 
density /Ar((-, Aj, Uj)j=i^,,,^p) with respect to /xai ® ■ ■ ■ ® /iAp- 

Step 2. We now assume that {Ai, . . . , Ap) are random and that their joint distribu- 
tions have a density with respect to the Lesbegue measure on "H^, where "Hat is the 
space of Hermitian matrices of size A^. In particular, the matrices have almost surely 
A^ distinct eigenvalues, see [IB]. The spectral decompositions of the previous step are 
measurable (see [TS], Section 5.3]) and, with the above notations {Aj,Uj) for eigenval- 
ues and eigenvectors of Aj, we can write the joint distribution of {Ai, . . . , Ap) in the form 
fifAr((Aj, t/j)j=i,...,p)dAiAjv(^i) • • • dAiAiv(Ap)d/iwAr(f^i) • • • dfiUj^{Up). The symbol /^a^ denotes 
the Lebesgue measure on = {{xi, . . . , xn)\xi < ■ ■ ■ < xn} and fiuN is the Haar measure 
on the set Un of unitary matrices of size A^. For any rii, . . . ,np > 0, one has 



E 



-Tr[A-o...oA^.] 



p 

N 



^^N^^N ii,...,jp = l 

X 5(Ar((Aj,[/j)j=i,...,p)d/iA^f(Ai) . . .dii^^{kp)<l^iuM{Ui] 



N 



Wp /, • • • K:Mi(h^^ A,),=i,...,p)d^A^(AO . . . d/iA^(Ap) 

ii,...,ip=l "^^N 



where is as in the previous step and 

/7v((Aj,i, Aj, Uj)j=i^,„^p)g]\r[{Aj, f/j)j=i,...,p)d/it/jy (f/i) 



For any ii, 



,ip = 1, . . . , N , we have 



,dfiUj^{Up). 



K!i ■ ■ ■ >^7^ihN{{Xj,i, Aj)j=i,...,p)d/iA^(Ai), . . . , d/iA^(Ap) 
A"^ . . . Xp^'h^j^' ' *'''(Ai, . . . , Ap)dAi, . . . , dAp, 



(«iv,ip)^^^^ . . . , Ap) is obtained by integrating with respect to the variables A^j, 
for ki ii, kp ip. We finally obtain 



where h'p^' 



E 



A^ 



Tr[A^^ o . . . o A"''] 



A"^ . . . , Xp^'h^^Xi, . . . , Ap)dAi, . . . , dAp, 



A. 
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where Tin = ^ Tin ip ^tv^' Hence, a family of random variables (ai, . . . , Cp) as in the 
proposition exists and its joint distribution has density Jin with respect to the Lebesgue 
measure on M^. 

Step 3. We now consider the general case. Let {Xi,...,Xp) be a family of indepen- 
dent random matrices, independent of [Ai, . . . , Ap), distributed according to the standard 
Gaussian measure on Hn with respect to the inner product {A,B) = NTt[AB]. By the 
regularizing process of convolution on Hermitian space, for any e > 0, the joint distribution 
of {AI, . . . , Ap) = {Ai + eXi, . . . ,Ap + eXp) has a density with respect to the Lebesgue 
measure. By the previous step, there exists a family of random variables (a^, . . . , ap such 
that E[(a^)"i . . . (ap"''] = Tr o ■ ■ ■ o (Ap"^]] for any ni, . . . , rip > 0. As e 

goes to zero, (af , . . . , ap converges in moments to a family of random variables as in the 
proposition. □ 
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